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Abstract 

In this paper, we discuss ahernative means of measuring the possible presence of local parity 
violation in relativistic heavy ion collisions. We focus on the phenomenon of charge separation and 
introduce the charged dipole vector Ql, which will measure the charge separation on an event-by- 
event basis. Using Monte Carlo events, we demonstrate the method and its discriminating power. 
In particular we show that such an analysis will reveal the strength of charge separation effect 
^ ■ and its azimuthal correlation with the reaction plane. We further show that our proposed method 

may be able to distinguish between the actual charge separation effect and effects due to certain 
^ [ two particle correlations. The connection to present measurements based on particle correlations 
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is discussed. 
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I. INTRODUCTION 

Topological objects in Quantum Chroniodynaniics(QCD) (and generally in non-Abelian 
gauge theories) have attracted persistent theoretical interests and are important in many 
aspects [l|. For example, instantons are known to be responsible for various properties of 
the QCD vacuum, such as spontaneous breaking of chiral symmetry and the Uyi{l) anomaly 
(see e.g. |2[|3|). Magnetic monopoles, on the other hand, are speculated to be present in the 
QCD vacuum in a Bose-condensed form which then enforce the color confinement, known 
as the dual superconductor scenario for QCD confinement which is strongly supported by 
evidences from lattice QCD (see e.g. |^|5|). Alternatively vortices are also believed to 
describe the chromo-electric flux configuration (i.e. flux tube) between a quark-anti-quark 
pair in the QCD vacuum which in turn gives rise to the confining linear potential (see e.g. 
reviews in [a] [61]). Some of these objects, such as monopoles [Tl] and flux tubes 8|, may 
also be important degrees of freedom in the hot and deconfined QCD matter close to the 
transition temperature T^, and may be responsible for the observed properties of the so 
called strongly coupled quark-gluon plasma |9|. Certain phenomenological consequences of 



such topological objects for relativistic heavy ion collisions have been studied in 10 |. 
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A particularly interesting suggestion by Kharzeev and collaborators |llHl6|on the di- 
rect manifestation of effects from topological objects is the possible occurrence of V- and 
CV-odd (local) domains due to the so-called sphaleron transitions in the hot dense QCD 
matter created in the relativistic heavy ion collisions. In particular, the so called Chiral 
Magnetic Effect (CME)[13|] predicts that in the presences of the strong external (electrody- 
namic) magnetic field at the early stage after a (non-central) collision sphaleron transitions 
induce a separation of charges along the direction of the magnetic field. Since the external 
magnetic field is perpendicular to the reaction plane defined by the impact parameter and 
the beam axis, one expects an out-of-plane charge separation. As a result positive charges 
are expected to preferentially go in one (out-of-plane) direction and negative charges in 
the opposite (out-of-plane) direction. In a given event, this charge separation results in a 
momentum space electric dipole which breaks parity. However, the dipole moment will be, 
with equal probability, parallel or anti-parallel to the magnetic field depending whether the 
Chiral Magnetic Effect is caused by a sphaleron or anti-sphaleron transition. Consequently, 
the expectation value of the dipole or, more precisely, of the scalar product of the dipole 



and the magnetic field, will vanish. 

For the aforementioned reasons the CME will not give rise to a non-vanishing expecta- 
tion value of a P-odd observable. However, the fact that parity is broken event-by-event 
should be reflected in the variance of a P-odd observable, which, however is a P-even ob- 
servable. Therefore, other, non-parity violating processes may contribute which need to be 
well understood. 

Very recently the STAR collaboration has announced the first experimental evidence of 
a possible local parity violation phenomenon at the Brookhaven's Relativistic Heavy Ion 
Collider(RHIC) [17]. STAR has measured the differences between the in-plane-projected 
and out-of-plane-projected 2-particle azimuthal correlations for both same-charge pairs and 



opposite-charge pairs, as proposed by Voloshin in 18|. The data indeed show very interesting 
charged-pair correlation patterns that depend on the charges (same/opposite), reaction plane 
(in-plane/out-of-plane), average pt of pairs, and colliding nuclei (Au/Cu). At first sight, some 
features of the data appear to be consistent with what has been expected from the local 
parity violation phenomenon. However, as shown in Ref.[l9|, certain aspects of the data 
appear to be puzzling at present: contrary to expectations from the Chiral Magnetic effect, 
the same-sign pairs show a negative in-plane instead of a positive out-of-plane correlation. 
Consequently, an interpretation of the data in terms of local-parity-violation would require 
a nearly exact cancellation for all centralities of correlations due to the Chiral Magnetic 
Effect and those due to ordinary correlations. In addition, further studies have proposed 



alternative contributions to the observed signals [20|, |2l| , and various possib 
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related to the local parity violation J22|] and the Chiral Magnetic Effect 

recently discussed. At present, therefore, the existence of local parity violation in heavy ion 

collisions has not yet been definitively established and further detailed and more differential 

analysis of the current observable as well as the development of alternative observables are 

necessary. 

The purpose of this paper is to propose and study in detail an alternative observable which 
specifically measures the charge separation predicted by the CME. To this end we generally 
investigate the charge separation effect as an intrinsic charge-dependent particle azimuthal 
distribution, and propose a way of measuring the magnitude of the charge separation effect 
and its orientation relative to the reaction plane. The orientation, if experimentally mea- 
sured, will be essential for an evaluation of the proposed local parity violation and CME 



interpretation. 

In the following, we will first introduce an intrinsic charge-dependent distribution rep- 
resenting the charge separation effect. We will then propose the charged dipole vector Ql 
analysis method and demonstrate its discriminating power by applying it to Monte Carlo 
events. We will also present a detailed study of the impact of certain two-particle correla- 
tions which may contribute to charge separation and cause ambiguity in the interpretation 
of the present STAR data. We will show that even with in the presence of these background 
correlations the measurement of Q1 together with the correlations proposed by STAR may 
be able to unravel the charge separation, its orientation and the various sources contributing 
to it. 

II. THE CHARGE SEPARATION EFFECT 

Let us begin by specifying what we mean by the charge separation effect in this pa- 
per. Consider the distribution of final state hadrons in the transverse momentum space 
as schematically shown in the FigJT] If the "center" of the positive charges happens to be 
different from that of the negative charges, then there is a separation between two types 
of charges which may be quantified by an "electric dipole moment" in the transverse mo- 
mentum space. Such a separation may arise either simply from statistical fiuctuations or 
may be due to specific dynamics, such as the Chiral Magnetic Effect. The later case is 

Py (arbitrary) 
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FIG. 1: A schematic demonstration of the proposed simultaneous analysis of Q1 and Q2 vectors 
in the same event. 



best discussed in the intrinsic frame, i.e. the frame defined by the direction of the reaction 
plane \l//j.p.. There, we can define an intrinsic charge-dependent single particle azimuthal 
distribution, which, besides a possible momentum-space electric dipole moment, also allows 
for the presence of elliptic flow. 

f^ (0, g) ex 1 + 2 V2 cos(20 - 2^r,p,) + 2qxdi cos(0 - ^c.s.) (1) 

Here q and (f) represent the charge and the azimuthal angle of a particle, respectively. The 
parameters V2 and di quantify the elliptic flow and the charge separation effect, while \l/c.s. 
specifles the azimuthal orientation of the electric-dipole (see Fig{T]) and \1'r.p the direction 
of the reaction plane. It is important to notice that an additional random variable x = il 
is introduced. This accounts for the fact that in a given event we may have sphaleron 
or ant i- sphaleron transitions resulting in charge separation parallel or anti-parallel to the 
magnetic fleld. Consequently the sampling over all events with a given reaction plane angle, 
"^R.p, corresponds to averaging the intrinsic distribution /^ over x, namely / =< /^ >x'^ 
1 + 2f2Cos(2(/) — 2\E'p.p.). Physically speaking this means that the charge separation (or 
electric dipole, being P-odd) flips and averages out to zero, thus causing the expectation 
value of any parity-odd operator to vanish. However, since (x^) = 1 the presence of an 
event-by-event electric dipole may be observable in the variance of a parity-odd operator. 

Provided an accurate identiflcation of the reaction plane is possible one may also redeflne 
the angles by setting \E'p.p. = and replacing "^c.s. by A\E'(7 5'. = \E'c.s. ~ ^k.p. in Eq.(IT]). For 
A^c.s. = 7r/2 the char ge sepa ration term takes the form of (isin(0) as expected from the 



Chiral Magnetic Effect [m]l8\. 



For measurements related to heavy ion collisions one may reasonably assume particle 
charges to be |g| = 1 which is the case for almost all charged particles e.g. charged pions 
and kaons, protons, etc. We also emphasize that the above distributions does not contain a 
directed flow term for either type of charges. Finally one may also consider a p^-differential 
formulation of the charge separation effect or charge separation effects associated with higher 
harmonics in the azimuthal angle 0. 

Let us next discuss how the above deflned charge- dependent intrinsic single-particle dis- 
tribution contributes to the charged particle correlations recently measured by the STAR 
collaboration in [17l]. Before doing so we note that there will likely be additional contribu- 
tions from two- and multi-particle correlations which we will not consider in this Section. 



The STAR collaboration has measured charge dependent two- and three-particle correlations 
[l7|. Specifically they considered 

(i) The two-particle correlation < cos(0j — (j)j) > for same-charge pairs (+ + / ) and 

opposite-charge pairs (H — ). The contribution to this correlator due to the charge-depended 
intrinsic single-particle distribution, Eq.(IT]) is: 

< cos(0i - (t>j) >++/— = dl (2) 

< cos(0i - (pj) >+_ = -dl (3) 

(ii) The three-particle correlation < cos(0i + 0j —20^) > for same-charge pairs {i,j = 

+ + / ) and opposite-charge pairs {i,j = H — ) with the third particle, denoted by index 

k, having any charge. The contribution to these correlators due to the distribution, Eq.([T]) 
turns out to be 

< COs(0i + (f)j - 20fe) >++/ — , k-any = ^2 dj COs(2 A^c.S.) (4) 

< COs(0i + 4>j - 20fc) >+_, k-any = -^2 dl COs(2 A^c.5.) (5) 

where "k-any" indicates that the charge of the 3-rd particle may assume any value/sign. 



The STAR collaboration has demonstrated [17|] that the above three particle corre- 
lator is dominated by the reaction plane dependent two-particles correlation function 
< cos{(f)i + (pj — 2\E'ij.p.) > and within errors they have found that 

< cos(0i + (j)j - 2(j)k) >= V2 < cos{(j)i + (j)j - 2^ij.p.) > (6) 

Based on the distribution Eq.([T]) we find the same relation between these correlation func- 
tions, since the reaction-plane dependent two-particle correlation is given by 

< cos{(j)i + (t)j - 2^^.p.) >++/— = dl cos(2 A^c.s.) (7) 

for same-charge pairs, and 

< cos(0i + 0, - 2^^.p.) >+„ = -dl cos(2 A^c.s.) (8) 

for opposite-charge pairs. 

The proposed Chiral Magnetic Effect corresponds to A'^c.s. = f • I^i this case 

< cos(0i + (f)j - 2*p.p) >++/— = -dl<0, (9) 

< cos(0i - (pj) >++/— = +dl> 0. (10) 



Thus, the two correlation functions are expected to be equal in magnitude with opposite sign. 
The STAR measurement, on the other hand finds them approximately equal in magnitude 
but with the same (negative) sign. This discrepancy, which is discussed in detail in |l9|, 
needs to be understood before any definitive conclusions about a possible charge separation 
effect can be drawn. One aspect is the effect of higher order correlations, which we have so 
far ignored. 

In the following Section we present an alternative observable which is sensitive to a 
potential charge separation effect. As we will discuss in detail, this observable will take into 
account additional correlations beyond those considered by STAR and thus may help to 
clarify the present situation. 

III. MEASUREMENT OF CHARGE SEPARATION BY Qf ANALYSIS 

In this Section we will present a method to directly measure the intrinsic charge-dependent 
distribution in Eq.(IT]). We propose to measure the charged dipole moment vector Q1 of the 
final-state hadron distribution in the transverse momentum space. The magnitude Q1 and 
azimuthal angle \1'^ of this vector can be determined in a given event by the following: 

Ql cos \E'5 = ^ gi cos 0i 

i 

Q^sin^^^^g^sinc^, (11) 

i 

where the summation is over all charged particles in the event with qi and 0j the electric 
charge^ and azimuthal angle of each particle. This method is in close analogy to the Qi 
and Q2 vector analysis used for directed and elliptic fiow (see e.g. |25|). In particular the 
elliptic fiow and the reaction plane orientation can be determined via the Q2 vector (with 
magnitude Q2 and azimuthal angle \1'2) as: 

Q2 cos 2\I'2 = y^ cos 20j 
i 



^ As already mentioned, one may reasonably assume \qi\ — 1. 



We emphasize that contrary to Q2 the charge dipole vector, Ql, incorporates the electric 
charge qi of the particles^. The angles \1/^ and \l/2 are determined event-by-event from a 
finite number of particles and are not to be confused with the idealized expectations "^c.s. 
and \E'r.p., although they should become identical in the limit of infinite multiplicity. 

With both the magnitude Ql and the azimuthal angle "$1 determined on an event-by-event 
basis, important information about the underlying physics can be revealed.^ In particular, 
the event-by-event distribution of the magnitude Ql may indicate if there is a physical charge 
separation effect beyond pure statistical fluctuations. Furthermore, the relative distribution 
of the angle "^l with respect to the reaction plane determined in the same event (e.g. via the 
Q2 analysis) can show to which extend the charge separation is correlated with the reaction 
plane with a specific angle A'^c.s = "^c.s — ^_r.p. (see Fig{T]). 

Next let us investigate the discriminating power of the proposed Ql analysis for a potential 
charge separation effect. To this end we employ Monte Carlo sampling to generate an 
ensemble of final state hadrons with equal numbers of 7r+ and tt^. For each particle we sample 
an azimuthal angle (pi of its transverse momentum pt. We do not sample the magnitude of 
the transverse momentum, pt, as the Ql analysis defined in Eq. Olip involves only the angle 

To mimic an event that may have elliptic flow and a possible charge separation effect, 
we sample the azimuthal angles for A^+ = 200 7r+ and A^-_ = 200 vr" according to the 
intrinsic distribution in Eq.([T]) ^ ^. The sampling parameter V2 specifies the magnitude of 
the elliptic flow and \E'_r.p. the orientation of the reaction plane, which is randomly chosen in 
each event. The parameters di and "^c.s. specify the magnitude of charge separation effect 
and the orientation of the charged dipole, respectively. Finally for each event we randomly 



^ Mathematical details regarding the observable Q^ can be found in the Appendix. 

^ After the completion of the present work, we became aware of the recent PHENIX efforts 27| to measure 



the distribution of the difference of ^sin((/)i) between plus and minus charges. This is related to the 
quantity Ql sin^'j;, which is part of the full information that can be extracted by the Ql analysis. 
** One may actually assign p^-dependent weight factor in the definition of Q? to maximize manifestation of 



the desired physical effect, as has been done in the V2 analysis (see e.g. [25|). 
^ To be more realistic one shall allow for fluctuations of the particle numbers, N+ and N-, within a selected 

multiplicity window. These additional sources of background statistical fluctuations bring in negligible 

broadening of the Ql magnitude distribution curves in Figl5] as we have verified numerically with up to 

±20% fluctuations independently for A^_|_ and 7V_. 
^ We emphasize again that in such sampling, no directed flow effect will be generated and only pure statistic 

fluctuations will contribute in the usual Qi analysis. 



pick the sign of the parameter x = ±1. 

In order to investigate different physical situations we consider a variety of cases as follows: 

Case-la with di = and V2 = 0, and Case-lb with di = but V2 7^ 0. In both cases the 
charge separation may arise only from statistical fluctuations with or without the presence 
of elliptic flow; 

Case-II with di ^ but f 2 = 0. Here we allow for an explicit charge separation in addition 
to statistical fluctuations but do not consider any explicit elliptic flow; 
Case-IIIa with di ^ and f 2 7^ but with the two orientations "^c.s. and '^r,p, randomly 
assigned in each event, representing a situation where both charge separation and elliptic 
flow are present but not correlated in the azimuth; 

Case-IIIb with di y^ and f 2 7^ and with the two orientations "^c.s. and "^r.p, always 
parallel in each event; 

Case-IIIc with di ^ and f 2 7^ and with the two orientations "^c.s. and "^r.p, always 
perpendicular in each event; 

Case-IIId with rfi 7^ and 1^2 7^ and with the two orientations ^c.s. and '^r,p, always j 
apart in each event. 

A comparison between the various cases will allow us to test the ability of the proposed 
analysis to discriminate between the various physical scenarios. 

For each of the above cases we sample 100 million events. For each event we perform the 
Q2 and Ql analyses and extract the values (^25^2 and Ql,'^l. The resulting distributions 
of the magnitude Ql and the angle "$1 for the charged dipole vector Ql will be presented 
and discussed in the following subsections. 

A. The Magnitude Distribution 

The Ql magnitude distributions (normalized to unity) for the various aforementioned 
cases are shown in Figj2l We have chosen the following values for the sampling parameters: 
For the cases with non-vanishing dipole moment, Case-II and Case-IIIa,b,c,d, we have used 
di = 0.05. For the cases with non-vanishing elliptic flow. Case-lb and Case-IIIa,b,c,d, we 
allowed for two values of elliptic flow parameter V2 = 0.1 and V2 = 0.05 in each case. 
Therefore, overall there are twelve curves in the plot: interestingly, they all fall into only two 
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FIG. 2: The Q^ magnitude distribution extracted from sampled events for the various cases. The 
curves for Case-la (red) and for Case-lb with two values of U2 (magenta) appear on top of each 
other at the left, while the curves for Case-II (green) and for Case-IIIa,b,c,d with two values of V2 
for each(blue) appear on top of each other at the right (see text for details). 

groups. One group includes all the curves corresponding to the cases with di = (left set of 
curves in Fig|2]) and the other corresponds to all cases with finite dipole moment, di = 0.05, 
irrespective of the presence and strength of elliptic flow. 

From these result we may conclude the following: (a) the Q1 magnitude distribution is 
very sensitive to the physical charge separation effect represented by di, but rather unaffected 
by the elliptic flow t>2; (b) the left curve(s) in Fig|2]represents the pure statistical fluctuation, 
and any deviation of the measured distribution from this curve would contain information 
about the physical charge separation effect. Experimentally the reference curve from sta- 
tistical fluctuation can be obtained by randomly re-assigning the charges of the particles in 
each event. We have verifled with our sampling, that the analysis after charge-reshuffling 
shifts all the curves to the statistical ones, denoted as Case la and lb in Figj2l 
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B. The Angular Distribution 

We next turn to the angular distribution. Since any experiment samples over all directions 
of the reaction plane the only experimental accessible information is the relative angle, 
A\E' = \E''j; — \E'2, between the direction of the charge dipole, \E'^, and that of the reaction 
plane, or more precisely, the elliptic flow, "^2- 

We focus on the four cases Case-IIIa,b,c,d, where both the elliptic flow and the charge 
separation effects are present. For studies in the present subsection, we set V2 = 0.1 and 
di = 0.05 in the sampling. For each sampled event we then extract the angles "$1 and \l/2 
from each event and calculate their difference A\E' = \E''[ — \I'2- Since the angle \E'2 (essentially 
the reaction plane orientation) is defined modulus vr, A\E' is equivalent to A^ ± tt and we 
may always transform A\E' to be in the interval A\l/ G [— 7r/2, vr/2)'^. Furthermore, what 
matters most for our discussion is the distinction between the in-plane and the out-of-plane 
orientations. Thus it is sufficient to know the absolute value of A\E', i.e. |A\E'| G [0,7r/2]. 

The resulting |A^| distributions for the Cases IIIa,b,c,d are shown in Fig|3I^left). They 
clearly exhibit distinctive patterns for the various cases. Out of the four cases, the Case-IIIa 
serves as a background since in this case the orientations of the elliptic flow and the dipole are 
uncorrelate(P . Any deviation from the Case-IIIa |A\E'| distribution indicates an azimuthal 
correlation between the dipole and the elliptic flow. This background may be removed by 
subtracting the results from Case-IIIa in our model. Such subtracted |A\E'| distributions for 
the Case-IIIb,c,d are shown in Figj3]^right) for comparison. 

The curve for Case-IIIb, where the orientations for the elliptic flow and charge separation 
effects are taken to be parallel, "^c.s. = "^r.p^ shows a maximum at |A\E'| = (in plane) and 
a rapid decrease toward |A\E'| = it/ 2 (out-of-plane). Thus the dipole Ql is predominantly 



^ In addition, the direction of the electric dipole moment is only known modulo tt due to the random factor 
X in Eq.(IT]) representing the equal likelihood of sphalerons and anti-sphalerons. 

® We note that even for the Case-IIIa where there is no a prior correlation in the sampling, the distribution 
is not entirely flat but slightly favors a smaller angular separation between "^1 and ^2- This is not 
surprising as both Ql and Q2 are determined from the same set of particles in a given event. Given the 
probability distribution, Eq.([T]), it is clear that for any angular sampling parameters ^c.5. and ^_r.p. 
the maximum of the probability distribution is located at certain angle ^max (and '^max + tt for the 
opposite charge) between the two. As a result both angles, \1/J and ^2, extracted from the sampled events 
tend to slightly align with the maximal angle "^max- In an actual data analysis, this background Case-IIIa 
angular distribution can be well approximated by the angular distribution from pure statistical background 
obtainable via the previously mentioned charge- reshuffling, as we have verified with our sampling. 

11 





-^ 






1.2 


^^Rand 




^. 


-*-Para 


^ A 


c 1 


\ 


-♦-Perp 


y^ 


2 


\ 


-^n|^ 


/ 


= 0.8 


\_^ 






"^ 


rr-^r~^. ^ 1 1 1 1 


« 0.6 
Q 




^'''*-t-t-j 


|0.4 


-^ 


«^ ^ 


K^^^^^ 


0.2 
n 









8 



4 

|A»P| 



8 



2 



c 
o 



i 0.25 

(0 

Q 
■o 

0) 

u 

(0 

S -0-25 
(0 




FIG. 3: (left) The distributions of the relative azimuthal angle |A^| between the charged dipole 
vector Q1 and the reaction plane for Case-IIIa(black star),b(red triangle), c(green diamond), d(blue 
box), (right) The subtracted distributions (i.e. the difference with respect to the Case-IIIa distribu- 
tion) for Case-IIIb(red triangle), Case-IIIc(green diamond), and Case-IIId(blue box) respectively. 
(See text for more details). 

oriented in-plane. The curve for Case-IIIc, where the orientations for the elliptic flow and 
charge separation effects are taken to be perpendicular, shows exactly the opposite. Finally 
the curve for Case-IIId, where the orientations for the elliptic flow and charge separation 
effects are sampled to be 7r/4 apart from each other, features a maximum around | A\E'| = 7r/4 
and a decrease toward both ends. The effect, however is not as prominent as in the previous 
cases and it may be more difficult to separate it from the background Case-IIIa. For these 
scenarios we have also varied values of both f 2 and di in our sampling and found the patterns 
to be qualitatively the same (see Section [IV] for details). 

From these studies, one can conclude that by simultaneously measuring the orientations 
of Ql and Q2 and examining their angular difference distribution, the azimuthal correlation 
between the two physical effects may be extracted. 
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C. An Alternative Measurement of Relative Orientation 

While measuring the relative orientation distribution as discussed in the previous sub- 
section contains the "full" information, it is experimentally rather demanding as it re- 
quires determination of reaction plane in each event. There is an alternative way to 
determine whether the charged dipole Ql is closer to the in-plane or out-of-plane direc- 
tion without the need of reaction plane. The idea is to measure the correlation function 
(cos(2A\E')) = (cos(2\E'c'.5 — 2\E'/j.p)),^. which in a sense represents the "elliptic anisotropy" 
of Ql. By using the definitions in Eqs.( iTT]|T2l) . one obtains 

' N,f, + 2{t,j}'i + {t,j}2 + {t,j; kY 
[Nch + {^,J}'i]■[Nch + {^,J}2] 
with the two- and three- particle correlations defined as 



(cos(2A^)) = ( ;!: . .V 's I. .1/2 > (13) 



{hjVi 


= Yl ^^^j ^^^^^^ ~ ^ 




i¥=j 


{^,J}2 


= y ]cos2{(j)i-(j)j) 




iH 


{^,j;kr 


= Yl ^*^j ^^^^^^ + 




i^j^k 



_ - 20fc) (14) 

(see the Appendix A for details). 

If the observable (cos(2A\E')) is measured to be unambiguously negative/positive then 
the charge separation is closer to out-of-plane/in-plane. If it is however consistent with zero 
within errors, then the situation is unclear. An ideal Chiral Magnetic Effect without any 
other effect and without statistical fluctuations would predict (cos(2A\E')) = — 1. As will be 
discussed in the next Section, however, statistical fluctuations and additional two-particle 
correlations may lead to sizable corrections. 

Let us finally point out that the correlation function, Eq.f lT3|) . in principle involves particle 
correlations to all orders, since it represents an average of a ratio. In addition, the correlators 
appearing in the denominator, in particular {i,j}2, is of the same magnitude as Nch, which 
makes an expansion of the denominator unreliable. Consequently, the correlation function, 
Eq.( !T3|) . has to be measured in its entirety and represents a different measurement than that 



^ A similar correlation between the directed flow and elliptic flow was proposed by Poskanzer and Voloshin 
in 12611 to determine their relative orientation. 
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of the individual terms, which has been already carried out by the STAR collaboration. A 
dedicated measurement is thus important and urgently called for. 

IV. DISENTANGLING CHARGE SEPARATION AND TWO-PARTICLE COR- 
RELATIONS 

In this Section we focus on the influence of certain two-particle correlations on both the 
observables proposed in the present paper and the correlations measured by STAR. As has 
been shown in |19|, the present STAR data indicates the existence of at least two types 
of correlations (beyond the possible correlations predicted by the Chiral Magnetic Effect): 
same-charge pair back-to-back correlations (mostly in-plane), and opposite-charge pair same- 
side correlations (about equally in-plane and out-of-plane) . It is therefore important to 
study the robustness of various observables against such "background correlations", and 
in particular to evaluate the prospects of disentangling a physical dipole charge separation 
from those background correlations. Physically speaking, there are a number of potential 
sources that may induce such two-particle correlations, such as clusters, rnomentum and/or 



local charge conservation [20|][2l|, or even chiral magnetic spiral effect [24 1. 

In order to study correlation effects we introduce correlations on top of our previous 
sampling based on single particle distribution in Eq.([T]). In particular we test the two types 
indicated by data, i.e.: (a) same-charge pair back-to-back correlation (SCBB), which we 
implement by randomly selecting a small fraction of same-charge pairs and sampling them 
according to |0i — 02| > 7r/2; (b) opposite-charge pair same-side correlation (OCSS), which 
we implement by randomly selecting a small fraction of opposite-charge pairs and sampling 
them according to |0i— 02| < 7r/2. At the same time we ensure the single particle distribution 
remains unchanged. ^° In the following we discuss various situations with different choices 
for the strength of the dipole charge separation and the correlations. For all studies in this 
Section we set f2 = 0.1 and generate 10 million events for each case. 



^'^ The precise procedure we used here is the following: when sampling each particle (say a plus charge) we 
trigger the correlation with probability 0.4 and once triggered we correlate it with the preceding plus/minus 
charge for SCBB/OCSS correlations by sampling its angle according to both the single particle distribution 
and the constraint from correlation. We have studied different implementations (which essentially imply 
different approximations to the desired correlations) and found that our conclusions are not changed. 
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A. Test-I 

In Test-I, we consider a setting with relatively "large" out-of-plane dipole and small 
correlations and study four cases: 

Case-IIIc (previously studied) with dipole di = 0.05 and |A\I'| = tt/2 in the sampling; 
Case-IIIc-(a) with the same dipole as in Case-IIIc, plus SCBB correlation at a fraction 
fscBB = 0.4% of all same-charge pairs ; 

Case-IIIc-(/3) with the same dipole as in Case-IIIc, plus OCSS correlation at a fraction 
focss = 0.4% of all opposite-charge pairs ; 

Case-IIIc-(7) with the same dipole as in Case-IIIc, plus both SCBB and OCSS correlations 
at the same fraction of fscBB = focss = 0.4% . 

The results for Q1 and |A\E'| distributions are shown in Fig JH and the various measured 
observables ^^ are reported in TableJH From these results one can see that: 

(1) Both types of correlations shift the Ql distribution towards smaller values, i.e. they 
suppress the charge separation; 

(2) The |A\E'| distribution, on the other hand, is quite robust in the present setting and 
and maintains a clear pattern characteristic of an out-of-plane charge separation in all four 
cases. However, the correlations lead to a slight flattening of the distributions. 

(3) Correspondingly, the observable (cos(2A\E')) remains negative for all cases, with its 
absolute value slightly reduced by both types of correlations; 

(4) The observable (cos(0i + (pj — 20^)) remains negative for same-charge pairs and positive 
for opposite-charge pairs; 

(5) The observable (cos((/)j — (/)j)) has completely different sign patterns in the four cases, 
and appears to be a sensitive diagnostic for different types of correlations; 

(6) The Case-IIIc- (7) (see last column in Table JT]), where the dipole coexists with both SCBB 
and OCSS correlations, is the only one which shows a sign patterns for all four measured 
charge correlations that is qualitatively similar to the STAR data.^^ 

These results can be partly understood as follows: the physical dipole separates positive 



^^ The definitions and interrelations of tliese observables can be found in the Appendix. 

^'^ One could imagine that tuning the strength of each of the three components may provide a reasonable 

fit to the presently available STAR data. Given the rather schematic correlations employed in this study, 

the value of such an exercise is, however, rather limited. 
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■B 0.04 





FIG. 4: The Qf (left) and |A^| (right) distributions for the four cases in Test-I. The triangle (black), 
diamond(red), box(green), and star(blue) symbols are for Case-IIIc, Illc-(a), IIIc-(/3), and IIIc-(7) 
respectively (see text for more details). 



TABLE I: Various obervables in Test-I (see 


text). 




(0) 


IIIc 


Illc-(a) 


IIIc-(/3) 


IIIc- (7) 


m?) 


799.8 


437.8 


437.6 


253.1 


{{Q2?) 


1998.2 


2024.8 


2027.9 


2047.0 


((AQ2).Q2) 


-13109.4 


-7345.3 


-7400.1 


-4314.1 


(cos(0, - (t)j))^+/__ 


0.0025 


-0.00078 


0.00197 


-0.00131 


{cos^cpi - 0j))+_ 


-0.0025 


-0.00125 


0.00149 


0.00053 


{cos{2cP, - 2</.,-))++/__ 


0.0100 


0.0102 


0.0102 


0.0103 


(cos(20i-2,/.j))+_ 


0.0100 


0.0102 


0.0102 


0.0103 


^ {cOs{(t>^ + (t>j - '^(t>k))++,__^k~any 


-0.0025 


-0.00167 


-0.00107 


-0.00109 


^(cos(</>i + <^,--2./.fc))_^_^,_^„^ 


0.0025 


0.00130 


0.00192 


0.00082 


(cos(2A^)) 


-0.3118 


-0.2873 


-0.2889 


-0.2636 



and negative charges in the out-of-plane direction, while both the SCBB and the OCSS 
correlations, either by separating same-charge pairs or by focusing opposite-charge pairs. 
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tend to reduce the out-of-plane charge separation. 



B. Test-II 



It is important to assess how the above features change with the magnitude of the built- 
in dipole. We thus also study a similar setting to Test-I but with a smaller out-of-plane 
dipole di = 0.02. Again, we have four cases: Case-IV with the smaller dipole, and Case- 
IV-(a),(/3),(7) with the smaller dipole plus the same correlations as in the Case-IIIc- 
(a),(/3),(7). The results are shown in Figj5] and TablellTI The qualitative conclusion is 
largely the same as in Test-I, except for a few differences: 

(1) Not surprisingly the angular distribution is much flatter due to a smaller dipole, and the 
(cos(2 A\E')) is getting much closer to zero though still negative; 

(2) Case-IV-(/3) is interesting, since one finds ^ (cos(0i + (j)j — 20^))^^, k-any ^° ^^ 

positive although there is still out-of-plane charge separation (as indicated by angular dis- 
tribution and negative (cos(2 A\E'))). 

C. Test-III 

In this test we completely turn off the dipole, i.e. di = 0, in order to expose the effects due 
to correlations. Here the Case-lb (previously studied) has elliptic flow but no dipole and 
represents charge separation from pure statistical fluctuation, while Case-Ib-(a),(/3),(7) 
have no dipole but the same correlations as in the Case-IIIc-(a),(/3),(7). The results are 
reported in Figj6]and Table JITTf ^. A comparison with Test-I, II shows the following: 

(1) The correlator (cos(0j — (pj)) for both same-charge and opposite charge pairs show sim- 
ilar patterns as before, which indicates that they are most sensitively dominated by the 
correlations; 

(2) Most surprisingly, we find the correlations to rotate the angular distribution in opposite 
direction to the previous cases; 

(3) Furthermore Case-Ib-(7) (see last column of Table JITTj) . even without a physical dipole, 
shows a sign pattern for all four measured charge correlations qualitatively similar to the 



^^ In Table Unl and ITVl certain extremely small values are given in the form of 6(10 "), which means they 
are consistent with zero provided our finite statistics and may exactly vanish at infinite statistics limit. 
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Case-IIIc-(7) and to the STAR data. 
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FIG. 5: The Q1 (left) and |A^[ (right) distributions for the four cases in Test-II. The trian- 
gle(black), diamond(red), box(green), and star(blue) symbols are for Case-IV, IV-(q), IV-(/3), and 
IV- (7) respectively (see text for more details). 

TABLE II: Various obervables in Test-II (see text). 



(^) 


IV 


IV- (a) 


IV-(/3) 


IV-(7) 


m?) 


464.0 


269.6 


269.8 


166.5 


{{Q2?) 


1996.5 


2001.4 


2001.7 


2004.5 


((AQ2).Q2> 


-423.8 


-704.1 


-716.5 


-719.3 


(cos(0i - 0j))++/__ 


0.0004 


-0.00184 


0.00091 


-0.00185 


(cos(0i - 0j))+_ 


-0.0004 


-0.00020 


0.00253 


0.00107 


{co^{2(t>i - 2(/>,))++/__ 


0.0100 


0.0100 


0.0100 


0.0101 


(cos(20, - 2</.,-))+_ 


0.0100 


0.0100 


0.0100 


0.0101 


^ (cos(0, + 0, - 2(t)k))^^i^_^k_,^y 


-0.0004 


-0.00056 


0.00003 


-0.00049 


4{cos(,/.i + 0,-2,^,))^_^,_^„^ 


0.0004 


0.00021 


0.00080 


0.00022 


(cos(2A^)) 


-0.0232 


-0.0375 


-0.0387 


-0.0527 



It appears that on top of the charge separations due to statistical fluctuations, the cor- 
relations studied here seem to actually suppress the in-plane separation while "enhancing" 
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FIG. 6: The Q\ (left) and |A^| (right) distributions for the four cases in Test-III. The trian- 
gle(black), diamond(red), box(green), and star (blue) symbols are for Case-lb, Ib-(a), Ib-(/3), and 
lb- (7) respectively (see text for more details). 

TABLE III: Various obervables in Test-III (see text). 



(^) 


lb 


lb- (a) 


Ib-(/3) 


Ib-(7) 


m?) 


400.1 


237.8 


237.7 


150.0 


{{Q2?) 


1996.7 


1996.2 


1995.2 


1996.1 


((AQ2).Q2> 


1998.3 


543.1 


548.6 


-56.34 


(cos((/>i - (t)j))^^i__ 


o(10-7) 


-0.00204 


0.00071 


-0.00196 


(cos(0i - 0j))+_ 


6(10-7) 


6(10-7) 


0.00118 


0.00053 


(cos(2</)i - 20,))++/— 


0.0100 


0.0100 


0.0100 


0.0100 


(cos(2</., - 2</.,-))+_ 


0.0100 


0.0100 


0.0100 


0.0100 


^ (COS(0, + </),■ - 'i<l)^))++/__^k-any 


a(io-7) 


-0.00036 


0.00024 


-0.00038 


^ (cos(0i + 0j - 2<i)k)) +__ ^k-any 


0(10-6) 


6(10-7) 


0.00059 


0.00011 


(cos(2A^)) 


0.0532 


0.0242 


0.0245 


-0.0038 



the out-of-plane separation to some extent. The reason for this behavior is the presence of 
eUiptic flow: while the two-particle correlations themselves are reaction-plane independent, 
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due to the elliptic flow there are (on average) more pairs close to the in-plane direction than 
pairs close to the out-of-plane direction. Therefore, the correlations are more effective in the 
in-plane direction. Since, as already discussed, the correlations considered here tend to sup- 
press the charge separation, the statistical charge separations are more suppressed in-plane 
than out-of-plane. As a result we observe an effective "enhancement" of out-of-plane charge 
separation. This subtle effect is proportional to the magnitude of the elliptic flow, V2, and 
the correlation strength. 

D. Test-IV 

The results in Test-III suggest an interesting question: could it be that in a certain 
parameter region, the same out-of-plane charge separation (e.g. the same |A\1'| distribution) 
can be produced either by a physical dipole (e.g. due to Chiral Magnetic Effect) or by any 
one of the SCBB or OCSS correlations? It turns out that this is possible. Here in Test-IV 
we provide an example: 

Case-V-(a;) with only a physical out-of-plane dipole di = 0.025; 

Case-V-(/3) with no physical dipole and only SCBB correlation at a fraction fscBB = 1%; 
Case-V-(7) with no physical dipole and only OCSS correlation at a fraction focss = 1%- 
In addition we also include the previously studied Case-lb with only statistical fluctuation 
as a reference for comparison. The results are presented in FigjT] and Table llVl A few very 
interesting features can be immediately seen: 

(1) All three cases V-(a)(/3)(7) show almost identical | A\l/| distribution and (cos(2 A\&)), i.e. 
three distinctive physical effects can lead to quantitatively similar "apparent" out-of-plane 
charge separation; 

(2) The physical dipole in V-(a), however, shifts the Q\ distribution toward larger values 
while the correlations in V-(/3)(7) shift the distribution toward smaller values, thus making 
the Q1 distribution also an important discriminating probe; 

(3) The observable — (cos(0j -|- (pj — 20^)) has almost the same value in V-(a;) with physical 
dipole and in V-(/3) with SCBB correlation and thus can not distinguish the two cases, but 
it differentiates the V-(7) with OCSS correlation; 

(4) The correlator (cos(0j — (f)j)) again appears to be a sensitive observable, giving different 
sign patterns for the three cases. 
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TABLE IV: Various obervables in Test-IV (see text). 



(^) 


lb 


V-(a) 


V-(/3) 


V-(7) 


mi?) 


400.1 


499.8 


89.4 


89.4 


{{Q2?) 


1996.7 


1995.9 


1996.2 


1995.8 


((AQ2).Q2> 


1998.3 


-1762.9 


-499.4 


-497.2 


(cos(0i - 0j))++/__ 


6(10-7) 


0.00062 


-0.00390 


0.00733 


{cos{(j)i - <l)j))^_ 


6(10-7) 


-0.00063 


6(10-7) 


0.01118 


(cos(20, - 2,/.,-))++/__ 


0.0100 


0.0100 


0.0100 


0.0100 


(cos(20, - 20,))+_ 


0.0100 


0.0100 


0.0100 


0.0100 


^ (C0S((/), + 0j - 2</-fc))++/— ,fc_anj, 


6(10-7) 


-0.00062 


-0.00059 


0.00343 


^(cos(,/.i + 0,-2</>fc))^_^,_^„^ 


6(10-6) 


0.00062 


6(10-8) 


0.00400 


(cos(2A*)) 


0.0532 


-0.06243 


-0.05996 


-0.05994 
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The most important lessons we learn from these tests and especially from Test-IV are: 
first, the proposed | A\l/| distribution and (cos(2 A\l/)) can readily reveal the geometry of any 
potential charge separation; 

second, an out-of-plane charge separation may arise due a physical out-of-plane dipole (e.g. 
with Chiral Magnetic Effect) or due to specific types of two-particle correlations (e.g. SCBB 
or OCSS); 

third, the unambiguous way to disentangle these different sources contributing to an out-of- 
plane charge separation, is to measure all the available observables including both the Ql 
magnitude and angular distribution and the STAR observables ^ (cos(0i + (pj — 20^)) and 
(cos(0j — (pj)) as well. 

V. SUMMARY 

To summarize, in this paper we have introduced and studied a method to experimentally 
extract a possible charge separation in the relativistic heavy ion collisions. Specifically we 
have proposed to measure the event-by-event distribution of the charged dipole vector, Q^. 
This analysis is able to determine the magnitude of the dipole as well as its azimuthal orien- 
tation with respect to the reaction plane. Using Monte Carlo events, we have investigated 
the sensitivity of this method for various scenarios, including the presence of additional two- 
particle correlations. We have shown that the combined information of the magnitude and 
direction of the dipole can distinguish between effects due to certain two particle correlations 
and those from e.g. the chiral magnetic effect. 

As in the case of the elliptic flow analysis, the proposed method may be refined by 
studying the Ql distribution differentially in transverse momentum or rapidity. Further 
improvements may introduce transverse momentum or rapidity dependent weight factors in 
the Ql extraction via Eq.f lTT]) . 

To conclude, we have demonstrated that our proposed method provides additional dis- 
crimination power over already existing measurements. This is essential in order to defini- 
tively determine the presence of local parity violation in heavy ion collisions. 
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Appendix A 

In this Appendix we show some details related to the charged dipole vector Ql analysis 
and discuss its relation to various charged particle correlations. 

The magnitude of Ql itself involves a square-root and, therefore, contains all multi- 
particle correlations, as is evident from a Taylor expansion: 



Ql = VNch+{t,jri = K 



ch 



'^ + 2i^('-''5-8ik«^''^''>'^-: 



(Al) 



with {i,j}l = '^ij.j <ii<ij cos{(j)i — (f>j). For example, the third term {{i,j}f) in the above 
expansion involves 2,3,4-particle correlations. This can be exphcitly evaluated to give: 

(Kj}i)'= iV,,(iV,,-l) 

+2 22 cos(0j - (j)j) + 22 cos 2(0i - (pj) 

+2 ^ Qiqj cos(0i + (pj - 20fc) 

+ ^ qiQjQkqi cos(0i + (j)j -(j)k- (pi) (A2) 

We notice that the 3-particle correlation term in the above is precisely the one proposed in 



18| and measured in [17l|. However, a full Ql analysis reveals more information than those 
expressed by the two and three particle distributions. First, the expression for the magni- 
tude of Ql, Eq. (lAip . involves higher order terms than Eq. (]A2p . Thus it involves correlations 
beyond three particles. Second, the determination of the magnitude and azimuthal orienta- 
tion of Ql necessarily involves all particles in an event and hence correlations involving more 
than three particles. Third, knowing the magnitude distribution of Ql allows to calculate 
the average of any moment < (Qi)" >, which contains n-particle correlations. As discussed 



m [28|j [29| in general the average of the n— th moment of an extensive observable like Ql will 
necessarily involve n— particle correlations. 

As an example, from the definition of Ql, Eq. (fTT]) . one can calculate the 2"*^ moment: 

(Qi)^ = y, Qi (cos^ (pi + sin^ (pi) + 2^ QiQj (cos (pi cos (pj + sin (pi sin (pj) 

i i^j 

= A^ch + ^ QiQj cos{(pi - (pj) (A3) 
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which is directly related to two particle azimuthal correlations for both same-charge and 
opposite-charge pairs. This relation is in analogy to the familiar relation for elliptic flow, 
i.e. 

{Q2? = N,h + Y, cos(20, - 20,) (A4) 

Next we derive Eq.( [T3|) . We first evaluate the quantity 

(AQ2) . Q2 = [{QD^ ■ cos(2A^)] ■ Q2 (A5) 

by using the definitions in Eq. (lllfl2p : 

[{Qlf ■ cos(2A^)] ■ Q2 = [{Qlf cos(2^i)] ■ [Q2 cos(2^2)] 

+ [(g5)^in(2*i)] ■[g2sin(2^2)] 
= A^ch + 2 ^ Qiqj cos(0i - (pj) 

+ ^ COS 2(0i - 0j) + ^ QiQj cos(0i + (pj - 20fc) (A6) 

Combining the above with Eq. flA3|A4l) . we obtain 

\{Q\f ■ cos(2AvI/)] . Q^ 



cos(2A^) 



Ar,;, + 2{^,jK + {z,j}2 + {^,j;A:}^ 



(AT) 



with the two- and three- particle correlations {«,j}i, {«,j}2, \^■,3\^Y defined in Eq.( !T4|) . 
The Eq. flT3|) is simply the event averaged version of the above expression. 
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